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The  Solar  Optical  Universal  Polaruneter  (SOU PI  on  Space! a b  2  collected  movies  of  solar  granulation  com¬ 
pletely  free  from  the  distortion  and  blurring  introduced  by  the  Earth's  atmosphere.  Individual  images  in  the 
movies  are  diffraction-limned  130  cm  aperturei  and  are  not  degraded  by  pointing  jitter  Ithe  pointing  stability 
was  07003  root  mean  squareL  The  movies  illustrate  that  Ihe  solar  3  minute  oscillation  has  a  major  role  in  the 
appearance  of  solar  granulation  and  that  exploding  granules  are  a  common  feature  of  the  granule  evolution. 
Using  three-dimensional  Fourier  filtering  techniques,  we  have  been  able  to  remove  the  oscillations  and 
demonstrate  that  they  dominate  the  temporal  autocorrelation  functions  of  the  granulation  pattern.  When  the 
oscillations  are  removed,  the  autocorrelation  lifetime  of  granulation  is  a  factor  of  2  greater  in  magnetic  field 
regions  than  in  field-free  quiet  Sun.  Using  a  technique  called  local  correlation  tracking,  we  have  been  able  to 
measure  horizontal  velocities  and  observe  flow  patterns  on  the  scale  of  meso-  and  supergranuiatton.  In  quiet 
region*  the  mean  flow  velocity  a  370  m  s'1,  while  in  magnetic  regions  it  is  about  125  m  s'1.  We  have  also 
found  that  the  root  mean  square  fluctuating  horizontal  velocity  field  in  quiet  Sun  increases  from  0.4$  to  1.4 
km  s*‘.  and  in  strong  magnetic  field  regions  it  increases  from  0.3  to  0.75  km  s'1  as  the  measuring  aperture 
decreases  from  4'  to  r.  Combining  the  results  from  temporal  and  spatial  autocorrelation  functions  with  the 
velocity  measurements,  we  conclude  that  the  decay  of  the  temporal  autocorrelation  function  is  due  as  much  to 
mo  non  and  distoruoa  of  granules  as  to  lifetimes  of  elements  in  the  pattern.  Attempts  to  follow  granules  using 
image-processing  techniques  yields  a  lifetime  for  individual  features  on  the  order  of  5  minutes.  The  size  dis¬ 
tribution.  lifetime,  and  proper  motions  are  also  derived  and  compared  with  the  measurements.  By  superposing 
the  location  of  exploding  granules  on  the  average  flow  maps,  we  find  (bat  they  appear  almost  exclusively  in 
the  center  of  mesognnulation  size  Sow  ceils.  The  density  of  exploding  granules  is  sufficient  for  their  expansion 
fronts  to  cover  a  typical  mesogranule  in  900  s.  Because  of  the  noouniformity  of  the  distribution  of  exploding 
granules,  the  evolution  of  the  granulation  pattern  in  mesogranule  cell  centers  differs  fundamentally  from  that 
in  boundaries  It  is  dear  from  (his  study  that  there  is  neither  a  typical  granule  nor  a  typical  granule  evolution. 
Even  after  the  solar  oscillations  are  removed,  a  granule's  evolution  is  dependent  on  local  magnetic  flux 
density,  its  position  with  respect  to  the  active  region  plage,  its  position  in  the  mesognnulation  pattern,  and 
evolution  of  granules  in  its  immediate  neighborhood. 
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ABSTRACT 

The  Solar  Optical  Universal  Polarimeter  (SOUP)  on  Spacelab  2  collected  movies  of  solar  granulation  com¬ 
pletely  free  from  the  distortion  and  blurring  introduced  by  the  Earth's  atmosphere.  Individual  images  in  the 
movies  are  diffraction-limited  (30  cm  aperture)  and  are  not  degraded  by  pointing  jitter  (the  pointing  stability 
was  07003  root  mean  square).  The  movies  illustrate  that  the  solar  S  minute  oscillation  has  a  major  role  in  the 
appearance  of  solar  granulation  and  that  exploding  granules  are  a  common  feature  of  the  granule  evolution. 
Using  three-dimensional  Fourier  filtering  techniques,  we  have  been  able  to  remove  the  oscillations  and 
demonstrate  that  they  dominate  the  temporal  autocorrelation  functions  of  the  granulation  pattern.  When  the 
oscillations  are  removed,  the  autocorrelation  lifetime  of  granulation  is  a  factor  of  2  greater  in  magnetic  field 
regions  than  in  field-free  quiet  Sun.  Using  a  technique  called  local  correlation  tracking,  we  have  been  able  to 
measure  horizontal  velocities  and  observe  flow  patterns  on  the  scale  of  meso-  and  supergranulation.  In  quiet 
regions  the  mean  flow  velocity  is  370  m  s'1,  while  in  magnetic  regions  it  is  about  123  m  s*‘.  We  have  also 
found  that  the  root  mean  square  fluctuating  horizontal  velocity  field  in  quiet  Sun  increases  from  0.4S  to  1.4 
km  s'1,  and  in  strong  magnetic  field  regions  it  increases  from  0.3  to  0.73  km  s'1  as  the  measuring  aperture 
decreases  from  4"  to  1”.  Combining  the  results  from  temporal  and  spatial  autocorrelation  functions  with  the 
velocity  measurements,  we  conclude  that  the  decay  of  the  temporal  autocorrelation  function  is  due  as  much  to 
motion  and  distortion  of  granules  as  to  lifetimes  of  elements  in  the  pattern.  Attempts  to  follow  granules  using 
image-processing  techniques  yields  a  lifetime  for  individual  features  on  the  order  of  3  minutes.  The  size  dis¬ 
tribution,  lifetime,  and  proper  motions  are  also  derived  and  compared  with  the  measurements.  By  superposing 
the  location  of  exploding  granules  on  the  average  flow  maps,  we  find  that  they  appear  almost  exclusively  in 
the  center  of  mesogranulation  size  flow  cells.  The  density  of  exploding  granules  is  sufficient  for  their  expansion 
fronts  to  cover  a  typical  mesogranule  in  900  s.  Because  of  the  nonuniformity  of  the  distribution  of  exploding 
granules,  the  evolution  of  the  granulation  pattern  in  mesogranule  cell  centers  differs  fundamentally  from  that 
in  boundaries.  It  is  clear  from  this  study  that  there  is  neither  a  typical  granule  nor  a  typical  granule  evolution. 
Even  after  the  solar  oscillations  are  removed,  a  granule’s  evolution  is  dependent  on  local  magnetic  flux 
density,  its  position  with  respect  to  the  active  region  plage,  its  position  h:  the  mesogranulation  pattern,  and 
evolution  of  granules  in  its  immediate  neighborhood. 

Subject  headings:  Sun:  atmospheric  motions  —  Sun:  granulation  —  Sun:  oscillations 


I.  INTRODUCTION 

In  recent  years  there  have  been  excellent  reviews  of  granula¬ 
tion  by  Wittmann  (1979),  Beckers  (1981),  and  Bray,  Lough- 
head,  and  Dun-ant  (1984).  However,  a  number  of  what  are 
thought  to  be  significant  properties  of  the  granulation,  for 
example,  lifetime,  size,  and  evolution  (Kitai  and  Kawaguchi 
1979;  Oda  1984;  Kawaguchi  1980;  Roudier  and  Muller  1987), 
are  still  under  active  discussion.  Although  the  majority  of 
astronomers  believe  that  granules  are  thermal  plumes  over¬ 
shooting  from  the  convection  zone  into  the  photosphere  (Bray, 
Loughhead,  and  Durrant  1984;  Nordlund  1985),  there  is  some 
evidence  that  they  are  turbulent  eddies  (Zahn  1987).  The  rela¬ 
tive  contribution  of  wave  and  convective  motion  to  the  veloc¬ 
ity  field  of  the  photosphere  is  still  controversial  (Beckers  1981 ; 

1  L.  Acton.  D.  Duncan.  M.  Finch,  Z.  Frank,  G.  Kelly,  R.  Lindgren.  M. 
Morrill  N.  Ogle  (deceased),  T.  Pope.  R.  Reeves.  R.  Rehse,  and  R.  Wallace. 
Lockheed  Palo  Alto  Research  Laboratory:!.  Harvey.  J.  Leibacber.  W.  Living¬ 
ston.  L.  November,  and  H.  Ramsey,  National  Solar  Observatory :  G.  Simon. 
Air  Force  Geophysics  Laboratory. 


Nordlund  1985;  Cox  1988).  At  present  it  is  not  possible  to 
answer  the  question,  “What  is  a  typical  granule?"  In  view  of 
the  current  status  it  is  not  surprising  that  there  exists  no  better 
discussion  of  granule  time  histories  than  that  given  by  Mehl- 
tretter  (1978),  which  was  based  on  images  from  the  Spektro- 
stratoskop  balloon  flight  of  1975. 

In  our  opinion  there  are  four  primary  reasons  for  the  current 
difficulties :  a  limited  number  of  good  time  sequences  of  granu¬ 
lation  images,  a  lack  of  good  repeatable  definitions  for  physical 
properties  in  an  evolving  flow  pattern,  limitations  on  storage 
capacity  and  computational  power  of  previous  computer 
image  processing  systems  for  handling  high-resolution  time 
sequences,  and  a  shortage  of  high-resolution  numerical  simula¬ 
tions  of  granulation. 

Although  quality  granulation  photographs  have  been 
obtained  for  nearly  a  century  (Janssen  1896),  until  recently 
there  were  very  few  time  sequences  of  high  uniform  quality. 
The  Stratoscope,  Spektrostratoskop.  and  Soviet  Stratoscopic 
Observatory  balloon  programs  represent  the  main  data  base 
for  the  reported  investigations  of  granule  evolution.  While 
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these  telescopes  produced  many  excellent  images  in  quiet  Sun, 
problems  with  focus  shift  and  pointing  stability  compromised 
their  time  sequences.  Because  of  the  variations  in  quality  and 
pointing,  studies  of  these  data  have  depended  largely  on  sets  of 
photographs  rather  than  movies.  Ground-based  observatories, 
in  particular  Pic  du  Midi,  have  produced  truly  outstanding 
images,  but  again  movies  are  rare.  In  the  past  year  both  the  Big 
Bear  Solar  Observatory  and  the  Swedish  Observatory  on  La 
Palma  have  produced  some  good  to  outstanding  movies,  and 
computer-intensive  analyses  based  on  the  techniques  devel¬ 
oped  for  handling  SOUP  movies  are  now  underway. 

,  The  lack  of  good  movies  has  aggravated  objectively  defining 
a  granule  and  granule  evolution.  The  traditional  statistical 
tools  for  characterizing  intensity  fluctuations  in  the  photo¬ 
sphere  sidestep  this  problem.  Unfortunately,  they  also  include 
other  physical  processes  and  instrumental  effects  unrelated  to 
granules.  Spatial  autocorrelation  functions  (ACFs)  and  their 
Fourier  transforms,  and  intensity  power  spectra,  are  all  com¬ 
promised  by  lack  of  knowledge  of  the  combined  modulation 
transfer  function  of  the  atmosphere,  telescope,  and  detector 
system  (Bahng  and  Schwarzschild  1962;  Schmidt,  Knolker, 
and  Schroter  1981 ;  Nordlund  1984).  The  SOUP  data  are  only 
partially  free  of  these  problems.  Temporal  ACFs  do  not 
separate  decorrelations  due  to  size  and  shape  changes  from 
proper  motions,  nor  do  they  distinguish  changes  due  to  the 
convective  flows  from  other  patterns,  such  as  oscillations  and 
waves. 

Another  class  of  techniques  use  pattern  recognition  algo¬ 
rithms  to  locate  granules,  identify  the  same  granule  at  different 
times,  and  derive  various  properties  of  the  set  of  granules  select¬ 
ed.  We  have  found  that  each  algorithm  has  its  own  peculiar 
sensitivities,  and  it  is  very  difficult  to  get  consistent  results  from 
different  algorithms.  For  example,  attempts  to  isolate  granules 
by  thresholding  techniques  are  sensitive  to  small-  and  large- 
scale  intensity  variations  introduced  by  wave  patterns.  Most 
granule  identification  processes  ultimately  are  based  on  subjec¬ 
tive  criteria  of  the  observer.  Because  these  criteria  are  not  uni¬ 
versal  or  even  well  defined,  it  is  very  difficult  to  compare 
observations  of  different  groups,  or  to  have  a  great  deal  of 
confidence  in  solar  cycle  variations  of  granulation  parameters, 
unless  a  single  group  has  reduced  all  the  observations. 

The  study  of  solar  images  using  modern  image  processing 
systems  requires  large  amounts  of  storage  capacity.  Good  gran¬ 
ulation  images  have  a  resolution  of  075  or  better.  If  three 
samples  are  taken  per  resolution  element,  each  10"  square 
region  of  the  Sun  implies  a  60  x  60  array,  or  3600  pixels.  A 
region  large  enough  to  include  a  supergranule,  40",  needs  a 
240  x  240  array.  For  convenience  256  x  256  pixel  arrays  are 
usually  studied,  and  this  requires  between  64  and  128  kilobytes 
of  computer  memory.  A  time  sequence  of  0.5  hr  with  10  s  time 
resolution  requires  between  12  and  24  megabytes  of  memory. 
Fourier  transform  operations  typically  require  4-8  byte  preci¬ 
sion  for  the  data.  Thus  performing  three-dimensional  Fourier 
transforms  involves  data  blocks  of  several  hundred  megabytes. 
Until  recently,  computing  capabilities  on  this  order  have  not 
been  generally  available  for  the  reduction  of  solar  data;  of 
course,  better  resolution  compounds  these  problems. 

Because  of  the  difficulties  in  defining  flow  patterns  and  the 
computational  effort  required  to  identify  any  order  in  a  movie, 
the  solar  observer  may  be  greatly  aided  by  good  theoretical 
guidance.  Unfortunately,  three-dimensional  modeling  is  expen¬ 
sive.  Currently  32  x  32  x  32  volumes  can  be  examined,  which 
are  really  not  large  enough  to  understand  the  granule-granule 


interactions  caused  by  the  larger  scale  flows  of  meso-  and 
supergranulation  and  by  wave  generation.  It  is  to  be  hoped 
that  the  64  x  64  x  64  volumes  now  being  used  for  these  calcu¬ 
lations  will  provide  more  guidance.  Nonetheless,  the  computa¬ 
tional  simulations  of  Dravins,  Lindegren,  and  Nordlund 
(1981),  Nordlund  (1985),  Steffen,  Ludwig,  and  Kruss  (1988), 
and  Hurlburt.  Toomre,  and  Massaguer  (1984,  1986)  clearly 
have  many  of  the  features  seen  in  solar  images. 

The  flight  of  SOUP  on  Spacelab  2  (Title  et  al.  1986)  pro¬ 
duced  nearly  6000  good  granulation  images,  but,  much  more 
important,  it  produced  several  movies  of  almost  30  minutes  in 
length.  These  movies  are  essentially  diffraction-limited  (30  cm 
aperture)  and  are  free  of  frame-to-frame  changes  in  image 
quality  and  pointing  jitter  (07003  root  mean  square  [rms]).  The 
movies  are  qualitatively  different  from  any  previous  ground- 
based  movies.  From  the  first  viewings,  it  was  clear  that  oscil¬ 
lations  played  a  major  role  in  the  appearance  of  the 
photospheric  surface,  that  exploding  granules  were  much  more 
common  than  have  previously  been  suggested  (Title  et  al. 
1987),  that  there  were  regions  of  ordered  flow  (November  et  al. 
1987),  and  that  magnetic  fields  affected  the  time  history  of 
granulation  (Simon  et  al.  1988). 

The  SOUP  team  was  prepared  to  handle  the  massive  data 
reduction  task  required  for  movies.  An  analog  laser  disk 
system  was  developed  for  interactive  display  of  movies  and 
movie  sequences.  A  large  amount  of  digital  magnetic  and  laser 
disk  capacity  was  available  for  the  immediate  and  final  data 
products.  Interactive  image-processing  tools  suitable  for  han¬ 
dling  large  sets  of  images  were  in  place,  and  a  small  staff  was 
available  for  the  digitization  of  the  images  and  the  manage¬ 
ment  of  tapes  and  disk  storage  space. 

Because  of  the  analysis  capacity  we  have  been  able  to  apply 
a  number  of  image-processing  algorithms  for  the  character¬ 
ization  of  granulation.  In  particular,  three-dimensional 
Fourier  filtering  to  remove  solar  oscillations  has  been  possible. 
In  §  III  we  show  that  solar  oscillations  dominate  the  temporal 
ACFs.  When  these  oscillations  are  removed,  it  is  clear  that 
quiet  and  magnetic  regions  differ  significantly  in  temporal 
history.  The  differences  between  quiet  and  magnetic  Sun  are 
also  visible  in  their  average  intensities  and  their  rms  intensity 
fluctuations  (§  IV)  and  quasi-steady  horizontal  velocities  (§  V). 
The  more  accurate  estimate  of  the  ACF  lifetime  obtained  from 
Fourier  filtered  images  coupled  with  horizontal  velocities 
obtained  by  local  correlation  tracking  lead  us  to  suggest  in  §  V 
that  the  temporal  change  in  the  ACF  is  dominated  by  displace¬ 
ments  rather  than  lifetime  of  the  individual  granules. 

The  traditional  technique  of  determining  lifetimes  is  to  iden¬ 
tify  granules  and  then  track  them  from  beginning  to  end.  In 
§  VI  this  has  been  done  automatically  with  two  different  identi¬ 
fication  algorithms.  These  automatic  procedures  yield  lifetimes 
that  are  much  shorter  than  might  be  expected  from  the  tempo¬ 
ral  ACF.  This  result  probably  is  caused  by  the  fact  that  gran¬ 
ules  have  internal  structure — multiple  bright  points  and  dark 
lanes — which  changes  during  their  evolution.  We  have  also 
measured  granule  densities,  sizes,  velocities,  and  spacings,  and 
the  distribution  functions  of  these  properties.  Our  results 
bracket  measurements  found  previously.  The  range  in  granule 
properties  found  by  different  algorithms  shows  that  granule 
characteristics  are  strongly  dependent  on  the  process  used  to 
define  them. 

In  §  VII  we  discuss  exploding  granules.  In  the  SOUP  data  at 
least  half  the  solar  surface  is  affected  by  the  circle  of  expansion 
of  an  exploding  granule.  Often,  a  granule  is  terminated  (e  g.. 
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•quashed)  by  a  neighboring  exploder.  Further,  these  exploders 
are  concentrated  in  meso-  and  supergranuie  interiors,  so  that 
(heir  effects  are  not  uniform.  Because  of  this  fundamental  non- 
uniform  distribution,  lifetimes  of  individual  granules  are  often 
determined  by  events  in  thrir  local  neighborhood,  which  in 
turn  depends  on  their  position  in  the  mesogranulation  pattern. 

Finally,  in  §  VIII,  we  summarize  our  results  and  stress  that 
there  is  not  a  single  “typical''  granulation  evolution  pattern, 
gather,  the  surface  of  the  Sun  has  brightness  fluctuations 
which  are  caused  by  a  combination  of  convection,  turbulence, 
waves,  and  magnetic  fields.  It  is  impossible  to  separate  these 
effects  in  a  single  image,  and  the  first  steps  in  effecting  the 
•operation  in  movie  data  are  reported  here.  The  evolutionary  ~ 
history  of  a  bright  feature  on  the  surface  of  the  Sun  depends  on 
where  it  is  born  in  a  preexisting  flow  field,  how  close  it  is  to  an 
exploding  granule,  and  how  much  magnetic  flux  is  in  its  vicin¬ 
ity. 


n.  OBSERVATIONS 

The  granulation  data  were  all  collected  on  Kodak  SO  2S3 
film  using  the  SOUP  instrument  during  the  flight  of  Spacelab  2 
(Title  er  al.,  1986).  The  original  images  are  140'’  x  25CT  and  are 
taken  in  a  thousand  angstrom  band  centered  at  6000  A.  The 
images  are  very  uniform  in  quality  and  are  distortion-free. 
Using  the  Spacelab  Instrument  Pointing  System  and  an  inter¬ 
nal  fast  guider,  image  stability  of  07003  rms  was  achieved 
during  quiescent  periods,  and  peak  disturbances  were  less  than 
0703. 

The  original  film  was  digitized  using  a  1024  x  1024  CCD 
camera  developed  for  the  High  Resolution  Solar  Observatory 
program  of  NASA.  The  film  was  mounted  in  a  high-precision 
animation  projector  with  a  double  integrating  sphere  light 
source  and  imaged  onto  the  CCD  using  a  high-quality  copying 
lens  mounted  on  a  x-y-z  stage.  By  moving  the  lens  and  stage, 
any  section  and  magnification  of  the  film  frame  could  be  select¬ 
ed  for  digitizing.  The  12  bit  digital  images  were  transferred 
directly  to  the  Lockheed  Palo  Alto  Research  Laboratory  VAX 
780  for  processing  and  analysis.  The  raw  digital  images  were 
corrected  for  both  flat  field  and  dark  current  Solar  rotation 
and  slow  drift  of  the  pointing  were  measured  by  cross- 
correlating  consecutive  images  in  a  movie.  These  effects  were 
then  removed  by  shifting  images  rigidly  using  bicubic  inter¬ 
polation.  A  more  detailed  discussion  of  the  analysis  process  is 
contained  in  Simon  et  al.  (1988). 

For  this  report  we  have  studied  four  digitized  subsections, 
two  quiet  Sun  regions  (34*  x  34*  and  40”  x  40”),  and  two 
regions  centered  on  the  pores  (34*  x  34*  and  60”  x  60”).  The 
pores  are  part  of  AR  4682  at  approximately  SIS,  W31 
In  —  0.75).  Because  of  limitations  of  the  processing  computer 
and  our  image  display  system,  all  digital  images  are  256  x  256 
pixels.  This  implies  OT134  and  07161  pixels  for  the  quiet  Sun 
regions,  and  07134  and  07234  pixels  for  the  pore  region,  respec¬ 
tively.  The  film  frames  were  separated  by  2  s  in  time  and  cover 
1650  s.  The  results  reported  here  are  all  from  orbit  110 
(19:10:35  to  19:38:07  GMT  on  1985  August  5).  For  most  of 
the  results  reported  here  we  have  used  data  with  10  or  20  s  time 
separations. 

In  addition,  we  have  used  magnetograms  collected  at  the  Big 
Bear  Solar  Observatory  (BBSO)  during  orbit  1 10.  The  BBSO 
magnetograms  were  supplied  in  digital  form  on  magnetic  tape. 
They  were  rescaled,  rotated,  and  rectified,  so  that  the  SOUP 
and  BBSO  images  could  be  accurately  coaligned. 


in.  TEMPORAL  autocorrelation  measurements 


We  define  the  correlation  lifetime  of  granulation  by  the  Me 
width  of  temporal  ACF,  where  the  temporal  ACF  of  the  inten¬ 
sity  is  defined  as 


ACF(r) 


<d/(x,  y,  t0)Sl(x,  y,  t0  +  r)> 
<d/2(x,  y,  t0 )> 


where  x  and  y  are  the  spatial  coordinates,  r0  is  the  time,  and  r  is 
the  time  separation.  The  brackets  indicate  an  average  over 
space  and  time,  and 

SI{x,y,t)  -  l(x,y,t)  -  </(x,  y,  r)> 


The  ACF  has  been  calculated  for  granulation  in  quiet  Sun 
regions  essentially  free  of  strong  magnetic  fields  and  in  the 
magnetic  field  regions  in  the  vicinity  of  the  active  region.  The 
locations  of  quiet  and  magnetic  areas  were  determined  with  the 
BBSO  magnetograms  by  carefully  registering  features  visible 
on  both  the  SOUP  and  the  BBSO  images,  namely,  the  sunspot 
and  pores.  After  superposition,  isogauss  contours  from  the 
magnetogram  were  overlaid  onto  the  SOUP  images  using  our 
image  display  system.  Small  image  areas  (24  x  24  pixels,  called 
“  boxes  ")  were  then  chosen  to  be  either  outside  (quiet  Sun)  or 
inside  (magnetic  areas)  of  the  contours.  The  sunspots  and  pores 
were  avoided.  Results  for  larger  areas  were  obtained  by  averag¬ 
ing  ACFs  from  several  boxes. 

Four  image  sequences  were  used,  comprising  three  different 
fields  of  view.  This  allowed  us  to  compare  our  results  from 
different  areas  of  the  Sun  at  different  image  scales,  and  also 
allowed  comparison  of  quiet  Sun  and  magnetic  areas  at  the 
same  image  scale.  Figure  1  (Plate  10)  shows  a  quiet  Sun  image 
(top)  and  a  pore  region  image  (bottom).  Indicated  on  both  are 
the  sample  boxes  used  to  calculate  the  ACFs  and  the  BBSO 
magnetogram  contours. 

The  ACFs  from  the  quiet  Sun  sequences  are  shown  in  Figure 
2.  These  results  are  the  average  of  16  boxes  for  the  07161 
pixel  * 1  sequence  (solid  line)  and  the  entire  field  of  view  for  the 
07134  pixel'1  sequence  (dashed  line).  The  24  x  24  pixel  boxes 
were  used  in  the  former  case  because  it  contains  regions  with 
significant  magnetic  flux,  which  were  avoided.  The  quiet  Sun 
correlation  lifetime  is  about  5  minutes.  Measurements  from  a 
number  of  other  studies  (Mehltretter  1978)  are  also  included  in 
Figure  2.  Previous  ACF  lifetime  measurements  of  granulation 
have  ranged  from  3  to  7  minutes,  with  an  average  of  about  6 
minutes  (Wittmann  1979).  As  can  be  seen  from  the  figure,  both 
the  shape  of  the  ACFs  and  the  determined  lifetime  of  the 
SOUP  data  are  in  agreement  with  the  previous  values. 

The  ACFs  from  the  individual  boxes  indicated  in  Figure  1, 
which  cover  only  36  arcsec2,  are  qualitatively  different  from 
their  sum.  Figure  3  (solid  ctaves)  shows  the  ACFs  generated 
from  four  of  these  boxes  in  regions  of  quiet  Sun.  These  ACFs 
do  not  drop  monotonically,  but  rather  exhibit  oscillations  with 
periods  of  3-6  minutes  (depending  on  the  sample).  This  sug¬ 
gests  that  solar  oscillations  are  affecting  the  correlation  life¬ 
time,  which  is  not  too  surprising,  since  they  are  strongly  visible 
in  the  granulation  movies. 

The  average  ACF  over  a  sufficiently  large  region  does  not 
oscillate  as  strongly  as  a  small  region.  This  is  exhibited  in 
Figure  4,  which  shows  a  three-dimensional  isometric  ACF 
plot:  correlation  coefficient  versus  tinte  lag  and  area.  To  create 
the  figure,  we  have  calculated  ACFs  from  a  nested  sequence  of 
square  regions  varying  from  12  to  400  arcsec2.  As  the  area  of 
the  Sun  used  for  the  ACF  measurement  increases,  the  oscil¬ 
lations  decrease. 
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Flo.  2. — Autocorrelation  functions  from  two  quiet  Sun  refloats  original  data.  Solid  Ime:  O' 1 61  pud'1  usage  sequence;  dashed  tine.  0"134  pixel 
sequence  Superposed  ate  autocorrelation  measurements  from  previously  published  results  for  comparison  (dors,  crosses,  open  circles ). 
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Fig.  J — Autocorrelation  functions  in  6*  *  6'  bows  of  quiet  Sun  from  an  original  image  sequence  Isolid  curvesi,  showing  the  effects  of  intensity  oscillations.  The 
dashed  curves  ate  from  the  same  four  hoses  after  three-dimensional  Fourier  Altering  of  the  image  sequence  to  remove  the  effects  of  oscillations 
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Fn.  4. — Autocorrelation  function  in  quiet  Sun  plotted  is  i  three-dimensional  surface,  illustrating  how  oscillations  tend  to  dampen  with  increasing  region  size. 
The  x-aais  is  time,  the  y-aiis  is  region  sire,  and  the  z-a*is  is  correlation  amplitude. 


To  remove  the  effect  of  oscillations,  we  have  applied  what  we 
call  a  subsonic  Fourier  filter.  A  raw  image  sequence,  which  can 
be  considered  to  be  a  single  three-dimensional  function  of 
intensity  (/*,  /,)  and  time  (r),  is  Fourier  transformed  into  a 
function  of  kx,  kt,  and  at.  The  subsonic  filter  is  defined  by  a 
cone, 

co 


in  (*,  a>)-space,  where  k  and  at  are  spatial  and  temporal  fre¬ 
quencies  and  v,  is  the  maximum  phase  velocity.  All  Fourier 
components  inside  the  cone  (i.e*  with  phase  velocities  less  than 
tji  are  retained,  while  all  those  outside  are  set  to  zero.  Then  a 
new  sequence  of  images  is  calculated  by  an  inverse  Fourier 
transform.  For  the  subsonic  filters,  values  of  e,  of  3,  4,  and  3 
fan  s' l,  well  below  the  7  km  s'1  sound  speed,  were  used.  For 
these  values  of  the  velocity  cone  is  totally  inside  the  region 
occupied  by  the  solar  f-  and  p-modes. 

A  subsonic  filtered  movie  exhibits  very  little  solar  oscillation, 
and  the  same  is  true  of  the  ACFs  formed  from  subsonic  filtered 
data.  The  dashed  curves  in  Figure  3  are  made  from  the  same 
areas  as  the  solid  curves  but  are  from  subsonic  filtered  (c,  «  3 
km  s'1)  rather  than  original  data.  Figure  3  shows  a  compari¬ 
son  of  the  ACFs  from  300  arcsec3  of  quiet  Sun  for  (a)  the 
original  data  and  (6)  the  subsonic  filtered  image  sequences.  The 
quiet  Sun  ACF  lifetime  of  the  subsonic  data  is  1.3  times  greater 
than  the  original — a  lifetime  of  410  versus  270  s. 


We  have  investigated  the  effects  of  our  choice  of  subsonic 
filter,  and  whether  an  abrupt  cutoff  had  any  important  conse¬ 
quences.  We  also  repeated  the  subsonic  filter  and  ACF  calcu¬ 
lations  for  values  of  p  of  4  and  3  km  s'1.  Regardless  of  the 
filter  selection,  the  ACF  lifetime  that  resulted  was  within  10% 
of  the  original  3  km  s' 1  filter  value. 

Figure  6  shows  how  the  ACF  of  the  original  data  changes 
with  spatial  resolution.  The  solid  curve  (a)  is  the  original  data 
ACF,  and  the  dashed  curve  (c)  the  ACF  calculated  from  data 
degraded  in  resolution  to  1”  by  smearing  the  original  frames. 
Smearing  was  performed  by  calculating  the  Fourier  transform 
of  each  image  in  a  sequence,  filtering  in  the  spatial  frequency 
domain  with  the  appropriate  circular  aperture  diffraction 
pattern,  and  then  performing  the  inverse  Fourier  transform. 
The  insignificant  difference  caused  by  loss  in  resolution  prob¬ 
ably  explains  the  consistency  of  ACF  results  obtained  by  differ¬ 
ent  observers  at  many  observing  sites  (see  Fig.  2).  Because 
neither  the  p-  nor  the  /-modes  have  significant  power  on  the 
scale  of  an  arcsecond,  blurring  the  data  should  have  little  effect 
if  the  ACFs  are  dominated  by  solar  oscillations. 

On  the  other  hand,  granulation  does  have  significant  power 
on  the  scale  of  an  arcsecond,  so  that  blurring  subsonic  data 
should  have  a  significant  effect.  Figure  6 d  also  shows  ACFs 
calculated  from  blurred  subsonic  filtered  data.  Here  the  life¬ 
time  of  the  degraded  data  (d;  dashed  curie)  is  doubled  com¬ 
pared  with  the  unblurred  subsonic  data  ( b :  dotted  curie). 
Taken  together,  Figures  1-6  demonstrate  that  the  solar  oscil¬ 
lations  make  a  dominant  contribution  to  the  temporal  ACF 
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Fw.  5.— Autocorrelation  functions  (a)  for  quiet  Sun  original  data  (Fig.  la),  <h)  for  the  same  dau  Fourier  filtered  to  remove  the  effects  of  solar  oscillations,  (c)  for 
magnetic  region  original  dau  (Fig.  1H  and  (d)  for  the  same  dau  after  Fourier  filtering. 
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FiO-  6. — Autocorrelation  functions  from  0*143  pixel * 1  quiet  Sun  sequence  (al  for  original  data,  (6)  for  same  dau  Fourier  filtered  to  remove  the  effects  of  solar 
oscillations,  (c)  for  same  dau  as  in  (a)  except  sneered  to  I  *  resolution,  and  id)  for  same  data  as  in  (h)  except  smeared  to  1 "  resolution. 
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and  that  the  characteristics  of  the  granules  become  apparent 
only  when  the  oscillations  are  removed.  It  also  suggests  that 
the  differences  between  quiet  and  magnetic  Sun  may  be 
masked  by  the  oscillations  present  in  unfiltered  data. 

There  have  long  been  arguments  for  the  idea  that  granula¬ 
tion  is  different  in  magnetic  regions.  This  was  evidenced  by  a 
reduction  in  contrast  in  magnetic  regions — “abnormal 
granulation”  (Dunn  and  Zirker  1973),  increases  in  granule 
number  density  in  magnetic  regions  (Macris  1931),  and  differ¬ 
ences  in  the  shape  of  the  line  bisector  in  photospheric  lines  in 
magnetic  as  compared  with  quiet  regions  (Bo net  et  al.  1984). 
However,  because  of  the  variations  of  granule  properties 
reported  by  various  observers,  this  evidence  has  not  been  con¬ 
vincing.  There  is  also  evidence  that  the  mean  bisector  shape 
varies  with  solar  cycle  (Livingston  1984),  and  that  perhaps  the 
granule  density  changes  with  solar  cycle  (Macris  and  Rosch 
1983). 

To  observe  the  effect  of  magnetic  fields  on  the  granulation 
pattern  ACF,  12  boxes  were  selected  inside  the  70  G  contour  of 
the  magnetograms  (see  Fig.  1,  bottom)  but  well  outside  of  pores. 
The  magnetic  flux  contained  within  these  12  regions  varied 
from  3  x  10”  to  8  x  10”  Mx.  Shown  in  Figure  7  are  the 
ACFs  formed  from  original  (solid  curves )  and  subsonic  filtered 
(dashed  curves)  data  for  four  of  these  boxes.  Figure  3  shows  a 
comparison  of  the  summed  (all  12  regions)  ACFs  from  original 
(c)  and  subsonic  filtered  (<f)  data  in  the  magnetic  areas.  The 
lifetime  in  the  magnetic  region  is  more  than  a  factor  of  2 
greater  in  subsonic  than  in  original  data— 890  versus  420  s. 


Bahng  and  Schwarzschild  (1962)  also  compared  the  ACFs  in 
quiet  and  magnetic  areas,  but  they  saw  no  difference. 

Figure  3  summarizes  the  effects  of  oscillations  and  magnetic 
fields  on  the  temporal  ACF.  The  original  data  (a)  in  quiet  Sun 
and  (c)  in  magnetic  regions  differ  somewhat,  while  the  subsonic 
filtered  data  shows  substantial  difference  between  (b)  the  quiet 
Sun  and  (d)  magnetic  areas.  The  results  are  independent  of  the 
details  of  the  Fourier  filter  used.  We  also  repeated  these  calcu¬ 
lations  on  a  second  digitized  region  surrounding  the  pores. 
This  region  was  34'  x  34'  in  size,  and  therefore  had  the  same 
pixel  size  as  the  quiet  Sun  region  above.  These  results  were  in 
excellent  agreement  with  those  shown  in  Figure  3. 

IV.  INTENSITY  FLUCTUATIONS 

The  mean  intensity  of  granules  in  the  SOUP  passband 
increases  in  magnetic  field  areas,  while  both  oscillations  and 
granule  scale  intensity  fluctuations  in  time  are  reduced.  The 
mean  intensity  and  the  rms  intensity  fluctuation  over  the  28 
minute  observing  interval  have  been  computed  in  the  pore  and 
quiet  Sun  sequences,  for  both  original  and  subsonic  filtered 
data.  Comparison  of  the  original  and  subsonic  filtered  mean 
intensity  images  shows  that  they  are  virtually  identical, 
because  the  oscillations  present  in  the  original  image  sequence 
average  out  during  the  28  minute  observing  period.  The  mean 
intensity  image  for  the  pore  region  clearly  shows  areas  where 
the  granules  are  brighter,  and  comparison  with  a  simultaneous 
BBSO  magnetogram  reveals  that  these  areas  have  magnetic 
fields  present. 


Fig  7. — Autocorrelation  functions  in  four  6"  *  6"  bones  of  magnetic  regions  from  original  l  solid  curves)  and  Fourier  filtered  [dashed  curves I  to  remote  the  effects 
of  solar  oscillations 
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Fn.  1— Plot  of  the  mean  value  of  intensity  in  the  photosphere  as  a  function  of  magnetic  field  strength.  The  intensities  are  from  an  average  over  83  images  of  the 
pore  region  spanning  28  minutes.  The  magnetic  field  strength  is  from  a  registered  BBSO  magnetogram. 


Figure  8  quantifies  this  result  by  showing  the  mean  intensity 
as  a  function  of  magnetic  field  strength.  Each  data  point  in 
Figure  8  is  calculated  by  fitting  a  Gaussian  to  the  intensity 
distribution  in  each  magnetic  bin,  the  smallest  of  which  con¬ 
tains  at  least  1500  image  pixels.  The  size  of  the  error  of  the 
mean  is  smaller  than  the  size  of  the  dots  in  most  cases.  The 
mean  intensity  at  zero  field  has  been  normalized  to  100.  Signifi¬ 
cantly,  as  soon  as  the  magnetic  field  strength  is  greater  than  the 
noise,  the  mean  intensity  of  granules  increases  significantly, 
and  remains  higher  until  the  field  is  strong  enough  for  pores  to 
form.  The  formation  of  pores  is  indicated  by  the  sharp  drop  in 
mean  intensity  at  -275G. 

Comparison  of  the  rms  intensity  fluctuations  from  the  orig¬ 
inal  and  subsonic  filtered  image  sequences  shows  that  the 
intensity  fluctuations  are  considerably  reduced  in  the  subsonic 
filtered  images.  Figure  9  shows  the  mean  normalized  rms  inten¬ 
sity  fluctuation  as  a  function  of  magnetic  field  strength  for 
original  data  (solid  curve)  and  subsonic  filtered  data  ( dashed 
curve).  Abo  shown  in  the  figure  are  the  values  of  the  rms  fluc¬ 
tuations  in  both  sets  of  quiet  Sun  data  (diamonds).  In  all  the 
movies  the  rms  intensity  fluctuation  peaks  at  zero  magnetic 
field,  and  then  declines  with  increasing  field  strength  until 
pores  form.  Note  that  the  point  at  which  Figure  9  shows  an 
upturn  is  the  same  magnetic  field  value  where  Figure  8  shows  a 
downturn.  Pores  have  significant  intensity  fluctuations  because 
any  change  in  a  pore's  size  or  shape  will  cause  an  intensity 
fluctuation  due  to  the  contrast  between  the  pore  itself  and  the 
surrounding  bright  granules.  These  variations  within  the  pore 
region  are  consistent  with  our  subjective  impression  of  the 
**  vigor  "  of  the  granulation  in  the  movies. 

The  BBSO  magnetogram  used  in  the  construction  of  Figures 
8  and  9  has  a  resolution  somewhat  worse  than  2".  This  limits 


the  accuracy  of  its  registration  to  the  SOUP  images,  and  our 
knowledge  of  the  spatial  distribution  and  strength  of  the  mag¬ 
netic  flux.  The  existence  of  the  relationships  shown  in  Figures  8 
and  9  therefore  suggests  that  the  presence  of  a  magnetic  flux 
tube  affects  an  area  of  granulation  considerably  larger  than  the 
diameter  of  the  flux  tube. 

The  diamonds  in  Figure  9  indicate  the  value  of  the  rms 
fluctuations  for  quiet  Sun.  The  average  fluctuation  for  quiet 
Sun  (subsonic  filtered)  is  about  10%  higher  than  the  corre¬ 
sponding  value  from  the  zero  field  regions  in  plage  regions 
surrounding  the  pores.  Therefore,  even  the  granulation  in  zero 
field  regions  in  plage  differs  from  that  for  field-free  quiet  Sun. 
Magnetic  effects  are  not  limited  to  filigree  in  intergranular 
lanes. 

V.  HORIZONTAL  FLOWS  VIA  LOCAL  CORRELATION  TRACKING 

It  is  apparent  from  the  movies  that  granules  both  move 
randomly  and  are  advected  by  meso-  and  supergranulation 
scale  flows.  We  have  measured  the  horizontal  flow  velocities  by 
a  technique  we  call  local  correlation  tracking  (November  et  al. 
1987).  Briefly,  a  coarse  rectangular  grid  of  points  is  chosen 
covering  the  image  field  of  view.  A  spatial  cross-correlation 
between  two  images  (typically  60  s  apart  in  time)  is  computed 
in  the  neighborhood  of  each  grid  point.  The  position  of 
maximum  cross-correlation  yields  a  local  displacement  vector 
for  each  grid  point.  We  interpret  these  displacements  as  hori¬ 
zontal  velocities.  The  neighborhood  of  each  grid  point  used  for 
the  comparison  is  defined  by  a  Gaussian  apodizing  function, 
and  the  full  width  at  half-maximum  (FWHM)  of  this  Gaussian 
represents  the  spatial  resolution  of  the  velocity  measurement 

We  have  produced  movies  of  the  flow  field  in  this  manner  b> 
applying  local  correlation  tracking  to  a  time  senes  of  image 
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MAGNETIC  FIELD  (GAUSS) 

Fn.  9. — Plot  of  the  mean  value  of  the  rms  intensity  fluctuation  in  the  photosphere  as  a  function  of  magnetic  field  strength,  for  original  data  (solid  cun*)  and  for 
Fourier  Altered  data  (dashed  cured.  As  in  Fig.  8,  the  rms  intensity  fluctuations  were  calculated  from  a  28  minute  long  SOUP  image  sequence  centered  on  the  pores. 
The  two  labeled  QS  indicate  the  mean  rms  fluctuation  at  zero  magnetic  field  for  a  SOUP  quiet  Sun  region. 


pairs.  Composite  movies  have  been  created  by  overlaying  the 
flow  field  on  the  original  and  subsonic  filtered  data.  These 
1  movies  clearly  show  that  granules  are  advected  by  the  flow 
I  field,  that  magnetic  regions  have  smaller  steady  flows  than 
quiet  areao,  and  that  rneso-  and  supergranular  scales  are  visible 
in  the  displacement  field  of  the  granules  (Simon  et  al.  1988). 
Averaging  the  velocity  over  the  duration  of  the  movie  shows 
these  larger  scales  clearly,  by  suppressing  flows  of  individual 
granules  and  effects  of  3  minute  oscillations.  Figures  10  and  1 1 
(Plates  1 1-12)  show  the  average  velocity  field  overlaid  on  (a)  an 
image  mid  (b)  the  divergence  of  the  average  velocity  field  for  the 
quiet  Sun  and  the  pore  region,  respectively.  The  two- 
dimensional  divergence  of  the  horizontal  flow  is  roughly  pro¬ 
portional  to  the  vertical  flow  (November  et  al.  1987). 

The  mean  flow  speed  (magnitude  of  velocity  vector)  is  370 
i  and  163  m  s'1  in  the  quiet  Sun  and  pore  region  sequences, 
respectively.  Figure  12  shows  the  histogram  of  the  speeds  in  the 
I  quiet  Sun  (left)  and  the  pore  region  (right).  Figure  13  shows  a 
scatter  diagram  of  the  flow  speed  versus  magnetic  field  strength 
in  the  pore  region.  The  velocities  for  Figures  10-13  were  calcu¬ 
lated  using  a  Gaussian  with  FWHM  of  4"  and  a  grid  spacing  of 
|  275.  This  mask  size  was  chosen  to  suppress  the  local  expansion 

and  collapse  of  individual  granules  but  still  show  meso-  and 
supergranule  scale  flows.  The  magnetic  field  strength  is  the 
average  from  an  8  x  8  pixel  area  centered  at  the  location  where 
each  horizontal  velocity  vector  was  calculated.  Also  marked  on 
Figure  13  is  the  mean  flow  speed  ( solid)  as  a  function  of  mag¬ 
netic  field  and  the  average  speed  in  quiet  Sun  (heavy  dot).  The 
mean  speed  is  strongly  dependent  upon  magnetic  field;  it 
j  varies  from  273  km  s'1  at  zero  field  to  100  m  s'1  for  field 

i  strengths  greater  than  173  G.  Significantly,  the  quiet  Sun 


average  flow  speed  is  33%  greater  than  the  flow  speed  in  zero 
field  areas  near  the  pores.  This  implies  that  magnetic  field  can 
affect  surface  flows  for  a  considerable  distance  outside  the  flux 
tubes  themselves.  Figures  12  and  13  show  clearly  that  there  is  a 
difference  between  quasi-steady  flows  in  quiet  and  magnetic 
regions  and  that,  like  the  intensity  fluctuations,  these  flows  are 
also  different  in  granulation  in  zero  or  low  field  areas  in  plage 
than  for  field-free  quiet  Sun. 

Besides  demonstrating  the  existence  of  quasi-steady  flow 
fields,  local  correlation  tracking  can  also  produce  estimates  of 
the  instantaneous  velocities  of  the  granules.  The  horizontal 
velocity  vector  at  a  single  grid  point  shows  fluctuations  with 
time  due  to  granulation,  oscillations,  waves,  and  random  noise 
(due  primarily  to  film  grain).  To  separate  the  noise  from  the 
solar  signals,  we  have  calculated  the  Fourier  power  spectrum 
of  the  x  and  y-components  of  the  velocity  in  all  the  points  in 
the  image.  From  the  average  power  spectrum,  the  white  noise 
level,  due  primarily  to  film  grain,  is  easily  identified;  the  rms 
fluctuating  solar  velocity  and  the  rms  noise  can  then  be  calcu¬ 
lated.  We  have  measured  the  velocities  using  a  set  of  Gaussian 
masks  of  FWHM  of  4',  2',  173  and  1".  The  smaller  masks  show 
more  of  the  random  motions  and  the  rates  of  expansion  and 
collapse  of  individual  granules,  in  addition  to  large-scale  flows 
and  oscillations.  Unfortunately,  the  rms  noise  increases  as  the 
width  of  the  mask  decreases,  so  that  it  is  not  possible  to 
measure  rms  velocities  at  the  resolution  limit  of  the  images. 
Shown  in  Figure  14  are  plots  of  rms  speed  as  a  function  of 
Gaussian  FWHM  for  quiet  Sun.  nonmagnetic  areas  of  the  pore 
region  (magnetogram  <  40  G),  and  medium  and  strong  mag¬ 
netic  areas  of  the  pore  region  (100-400  G  and  >400  G. 
respectively). 
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Fig.  12. — Histogram*  showing  the  distribution  of  measured  flow  speeds  in  the  quiet  Sun  region  (left)  and  pore  region  (right) 


The  speeds  plotted  in  Figure  14  represent  the  fluctuating 
flow  at  a  single  point,  as  measured  in  different  regions  and  with 
different  spatial  resolutions,  the  sharp  increases  at  finer 
resolution  must  be  due  to  granulation,  since  most  power  in  the 
5  minute  oscillations  is  at  larger  spatial  scales.  This  is  borne 
out  by  movies  of  the  flow  vectors  overlaid  on  the  images,  in 
which  exploding  granules  are  especially  obvious.  In  strongly 
magnetic  areas,  the  velocities  are  systematically  45%  less  than 
quiet  Sun  values.  Thus,  magnetic  fields  strongly  suppress  the 
vigor  of  the  granulation  on  1*  scales.  Nonmagnetic  locations 


within  the  plage  are  essentially  the  same  as  quiet  Sun  far  from 
strong  fields,  which  implies  that  the  magnetic  fields  do  not 
affect  the  vigor  of  granular  flows  more  than  a  few  arcseconds 
beyond  the  boundaries  of  the  flux  tubes  themselves. 

From  Figure  14,  the  rms  speed  in  quiet  Sun  is  at  least  1.4  km 
s  "  *.  This  velocity  is  quite  high  compared  with  previous  mea¬ 
surements  of  the  vertical  component  of  photospheric  velocities  i 
(Beckers  1981).  It  is  possible  that  the  local  correlation  tracking  ! 
technique  is  measuring  phase  velocities  as  well  as  mass  motion. 
However,  the  convective  component  of  the  flow  is  driven  by 


MAGNETIC  FIELD  STRENGTH  (GAUSS) 

Fio.  13. — Scatter  plot  showing  the  flow  speed  versus  local  magnetic  field  strength  in  the  pore  region  (07134  pixel' 1  sequence).  Also  indicated  is  the  mean  flow 
speed  (solid  line)  as  a  function  of  field  strength  and  the  mean  value  from  the  quiet  Sun  region  (dot  labeled  QS  at  0  Gl. 
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Fig.  14. — Plot  showing  rms  velocity  vs  FWHM  of  the  Gaussian  used  in  the  local  correlation  technique.  Data  from  different  types  of  regions  are  designated  by 
open  circles  (quiet  Sun),  plus  signs  (plage),  crosses  (medium  magnetic  plage),  and  asterisks  (high-strength  magnetic  plage). 


thermal  buoyancy,  and  radiative  losses  should  cause  the 
brightness  front  to  lag  the  gas  motion.  Numerical  simulations 
by  A.  Nordlund  (1988,  private  communication)  verify  that 
brightness  fronts  lag  flows.  Thus,  our  velocities  may  even 
underestimate  the  true  values.  Details  of  the  relations  of 
brightness  and  mass  flows  will  be  discussed  in  a  separate  paper. 
Nonetheless,  the  velocities  measured  here  are  the  correct  veloc¬ 
ities  for  understanding  how  the  temporal  ACFs  evolve. 

The  temporal  ACFs  compare  sets  of  spatially  aligned  data 
sets  as  a  function  of  time  and  yield  information  related  to 
feature  lifetime.  Spatial  ACFs  compare  the  same  image  with 
different  spatial  offsets  and  yield  information  on  the  spatial 
scale  of  the  image.  The  spatial  ACF  is  defined  as 


ACF/A) 


(6l( x,  y,  t)Sl{x  +  A,  y,  t)> 
<3/2(jc,  y,  t)> 


where  A  is  the  spatial  offset.  The  spatial  ACF  can  be  used  to 
estimate  the  effect  of  flows  on  the  ACF  lifetimes  because  the 
displacements  can  be  considered  to  be  the  results  of  a  flow.  A 
velocity  v  which  lasts  for  a  time  t  will  cause  a  displacement 
A  —  vi.  We  define  a  displacement  lifetime,  Tt,  as  the  time 
required  for  the  rms  flow  pm  to  cause  the  spatial  ACF  to  drop 
to  1/e: 


ACF/o^/T,)  -  1/e. 

Figure  15  shows  (a)  the  spatial  ACF  for  subsonic  filtered  quiet 
Sun  and  ( b )  the  spatial  ACF  for  magnetic  regions.  From  the 
figure  the  1/e  width  is  450  km  and  500  km  for  quiet  and  mag¬ 
netic  regions,  respectively,  which  is  consistent  with  previous 
measurements  (Wittmann  1979).  On  average,  structures  in  the 
quiet  and  magnetic  regions  are  about  the  same  size  in  the 
SOUP  images.  While  we  have  shown  that  granules  in  the  mag¬ 


netic  region  are  different  from  quiet  g,anules,  the  spatial  ACF 
is  not  a  good  discriminator  of  that  difference.  Nevertheless,  the 
eye  can  see  a  subtly  different  texture  in  the  magnetic  parts  of 
the  pore  region  in  single  images. 

Either  random  or  locally  ordered  motions  will  cause  the 
lifetime  measured  from  the  temporal  ACF  to  underestimate  the 
lifetime  of  the  individual  features  in  the  granulation  firld. 
Using  a  velocity  of  1.4  km  $  " 1  for  quiet  Sun,  the  displacement 
lifetime  T4  is  320  s,  compared  with  an  ACF  lifetime  of  270  s.  In 
the  medium  magnetic  plage  (see  Fig.  14),  we  use  1.0  km  s~  \ 
which  yields  a  T4  of  500  s,  compared  with  420  s  from  the  ACF. 
Thus,  in  both  cases,  the  flows  are  sufficient  to  cause  much  of 
the  decay  in  the  temporal  ACF.  Note  that  this  flow  includes 
not  only  the  displacements  of  entire  granules  but  also  some  of 
the  distortions  and  explosions  which  are  seen  by  the  local 
correlation  tracking  with  a  1*  FWHM  Gaussian.  The  temporal 
ACF  does  not  simply  reflect  the  lifetime  of  isolated  intensity 
features  in  granulation:  motion  and  distortion  of  the  features 
have  a  large  effect. 


VI.  GRANULE  TRACKING 

Because  of  random  and  ordered  motion,  the  temporal  ACF 
yields  a  poor  estimate  of  the  lifetime  of  individual  features  in 
the  granulation  pattern.  The  direct  method  for  determining 
lifetime  is  to  identify  and  then  track  individual  features  in  the 
pattern.  The  standard  approach  has  been  to  find  a  granule  on  a 
master  frame  and  then  to  track  it  as  long  as  possible  on  follow¬ 
ing  and  preceding  images  (Mehltretter  1978).  The  weakness  of 
all  human  identifications  is  the  biased  selection  of  the  tracked 
elements  and  the  determination  of  the  beginning  and  end 
points.  This  is  especially  true  because  the  majority  of  granules 
are  bom  from  an  existing  fragment  (Dialetis  and  Maoris  1984; 
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Mchltrettcr  1978).  Since  our  images  are  free  of  distortions 
caused  by  atmospheric  turbulence,  we  have  tried  to  develop 
unbiased  computer  procedures  for  identifying  and  then  follow¬ 
ing  granules. 

lie  first  step  toward  tracking  granules  is  to  identify  all  of 
them.  We  have  tried  two  different  algorithms,  which  we  call 
“center  finding"  and  “lane  finding."  The  first  is  based  on 
finding  locally  brightest  points  and  identifying  them  with 
granule  centers,  and  the  other  on  finding  the  intergranular 
lanes  which  surround  each  granule.  In  other  words,  center 
finding  finds  granule  centers  directly  by  identification,  and  lane 
finding  finds  their  boundaries  and  infers  a  center  in  the  middle 
of  an  enclosed  area.  Both  these  families  of  methods  suffer  from 
the  problem  that  a  “granule"  may  have  an  internal  structure 
with  multiple  bright  points  or  multiple  dark  lanes  during  its 
lifetime. 

All  center-finding  algorithms  are  based  on  exploring  a  local 
neighborhood  and  finding  the  bright  points  in  that  area.  The 
simplest  of  these  methods  is  based  on  thresholding— defining 
granules  as  the  connected  areas  above  some  intensity  thresh¬ 
old.  Pure  thresholding  methods  are  very  sensitive  to  intensity 
variations  over  the  frame  caused  by  seeing  and  solar  oscil¬ 
lations.  To  make  them  relatively  stable,  a  spatial  filter  must  be 
applied  to  the  images  to  remove  very  high  and  very  low  spatial 
frequencies  (Roudier  and  Muller  1987).  We  have  experimented 
with  thresholding,  but  we  feel  that  it  is  too  sensitive  to  low 
spatial  frequency  intensity  variations,  and  to  the  detailed  form 
of  the  spatial  filter  used. 

Instead,  our  center-finding  algorithm  searches  for  local 
maxima  by  comparing  the  brightness  of  every  point  in  the 
image  with  its  neighbors.  The  test  depends  upon  two  param¬ 
eters.  which  are  the  value  of  the  intensity  delta  by  which  a 
center  must  exceed  its  neighbors  and  the  size  of  the  neighbor¬ 
hood  examined.  Shown  in  Figure  16  (Plate  13)  is  (a)  a  typical 
granulation  frame  and  (b)  the  same  frame  with  the  centers 
marked.  The  intensity  delta  has  been  adjusted  to  give  a 


satisfying  result  for  (b).  The  number  of  centers  found  is  reason¬ 
ably  insensitive  to  the  size  of  the  local  neighborhood  examined 
in  a  range  from  074  to  079  (3-7  pixel  width).  The  average 
number  of  granules  identified  by  this  algorithm  is  500  in  a  1600 
arcsec2  area  of  quiet  Sun,  or  31  granules  per  100  arcsec2,  about 
the  mean  of  values  found  by  previous  observers  (Wittmann 
1979).  Based  on  this  statistic,  the  center-finding  algorithm  is 
finding  most  of  the  larger  granules,  but  not  the  fine  structures 
within  them,  and  not  the  smallest  features. 

The  lane-finding  algorithm  is  much  more  sophisticated.  It 
isolates  local  bright  areas  by  surrounding  them  with  a  bound¬ 
ary  lane.  The  granule  location  is  then  defined  as  the  center  of 
gravity  of  each  individual  hole  in  the  lane  net.  Figure  16c 
shows  the  quiet  Sun  frame  overlaid  by  the  lane  net  with  the 
granule  locations  marked.  About  1000  granules  were  found  by 
the  lane-finding  procedure,  or  63  granules  per  100  arcsec2, 
which  is  in  the  range  found  by  Namba  and  Diemel  (1969)  and 
Karpinsky  (1980). 

A  comparison  of  the  center-  and  lane-finding  results  is 
shown  in  Figure  16d,  where  the  granulation  mesh  is  shown 
filled  in  in  white  wherever  the  center-finding  algorithm  failed  to 
identify  a  granule  found  by  the  lane-finding  algorithm.  (The 
instances  where  a  single  granule  found  by  the  lane  finder  was 
identified  as  two  or  more  granules  by  the  center  finder  is  very 
rare).  By  inspection,  we  see  that  the  lane  finder  located  small 
features  missed  by  the  center  finder,  and  also  subdivided  some 
of  the  larger  granules.  The  granule  density  in  the  white  area  of 
Figure  16d  is  140  per  100  arcsec2.  With  the  assumption  that 
granules  are  in  a  roughly  hexagonal  close-packed  pattern,  the 
granule  center-to-center  distance  is  179,  1735,  and  079  for  den¬ 
sities  of  31,  63,  and  140  granules  per  100  arcsec2.  Again  these 
numbers  bracket  the  range  found  by  previous  observers 
(Wittmann  1979). 

The  results  of  the  lane-finding  algorithm  can  be  used  to 
measure  the  area  of  each  granule,  and  this  allows  us  to  gener¬ 
ate  a  number  of  statistical  properties  similar  to  those  generated 


No.  1.  1989 


STATISTICAL  PROPERTIES  OF  SOLAR  GRANULATION 


487 


by  Roudkr  and  Muller  (1987).  Our  areas  are  that  of  the  mesh 
that  isolates  a  granule.  We  will  define  the  “  granule  "  size  to  be 
a  circle  that  has  two-thirds  the  area  of  the  mesh.  We  do  this  so 
that  the  area  of  the  granules  and  that  of  the  intergranular  lanes 
are  approximately  equal.  The  exact  ratio  of  granule  to  lane 
area  has  little  effect  on  the  character  of  the  statistical  relations. 

Figure  1 7  shows  a  histogram  of  granule  number  versus  size, 
for  both  lane-found  ( solid  line)  and  center-found  ( dashed  line ) 
granules.  The  number  of  lane-found  granules  increases  mono- 
tonicaily  toward  smaller  size  until  a  cutoff  is  reached  due  to  the 
resolution  limit  of  the  SOUP  telescope.  This  means  that  SOUP 
cannot  resolve  all  of  the  granules  present,  a  fact  confirmed  by 
higher  resolution  ground-based  observations  that  reveal  struc¬ 
tures  smaller  than  the  SOUP  resolution  limit.  Figure  18  is  a 
histogram  of  fractional  area  occupied  by  granules  of  a  given 
size  versus  size.  Both  curves  show  a  broad  peak  between  078 
and  175.  in  agreement  with  Roudier  and  Muller  (1987).  Figure 
19  contains  the  distribution  function  of  center-to-center  dis¬ 
tances.  The  peak  of  this  distribution  is  determined  mainly  by 
the  observed  number  of  granules  per  unit  area  and  is  nearly 
independent  of  the  identification  algorithm  used.  The  area  dis¬ 
tribution  (Fig.  18)  is  a  reasonable  representation  of  the  length 
scale  of  granulation.  Even  if  there  may  be  a  large  number  of 
unresolved  small  granules,  their  contribution  to  the  total  area 
is  small,  so  that  it  is  fair  to  say  that  there  is  a  characteristic 
granule  size  in  the  vicinity  of  172—174. 

However,  we  see  no  evidence  for  the  abrupt  change  in 
granule  properties  at  173  reported  by  Roudier  and  Muller 
(1987).  Since  the  change  occurs  in  the  largest  granules,  the 
resolution  difference  should  not  be  critical.  From  inspection  of 
their  filtered  images,  we  believe  that  the  change  is  due  to  the 
effect  of  their  high-pass  filter  and  threshold  processing  on  the 
largest  granules.  Large,  relatively  smooth  intensity  features 
develop  long,  sinuous  perimeters  which  are  artifacts  of  the 
identification  algorithm. 


Given  a  set  of  images  with  the  granule  centers  marked,  it 
would  seem  straightforward  to  identify  the  center  of  the  same 
granule  on  sequential  movie  frames.  In  fact,  because  of  the 
tendency  for  granules  to  develop  multiple  centers  or  multiple 
lanes,  the  objective  identification  of  the  beginning  and  end  of  a 
granule  is  difficult  both  conceptually  and  computationally.  For 
our  algorithm,  identification  occurs  if  granule  centers  on  suc¬ 
cessive  images  are  separated  by  less  than  a  tolerance  A r.  By 
examining  the  entire  movie  starting  on  the  first  frame  in  which 
identification  is  possible  and  then  continuing  until  identifica¬ 
tion  is  no  longer  possible,  we  obtain  sets  of  x,  y,  and  t  coordi¬ 
nates  of  associated  granule  centers.  Each  association  forms  a 
trajectory  in  three-dimensional  spacetime,  which  we  call  a 
“  string,"  from  which  we  can  measure  the  velocity  and  lifetime. 
Because  the  string  algorithm  looks  at  each  frame  in  succession, 
it  will  not,  in  general,  properly  connect  a  granule  that  frag¬ 
ments  and  then  re-forms. 

Figure  20  contains  six  plots  of  the  number  of  granules  with 
lifetime  t  as  a  function  of  t.  It  illustrates  both  center-found  and 
lane-found  granules  in  both  original  and  subsonic  filtered  data, 
for  A r  of  075  and  1”.  All  of  the  plots  are  well  fitted  by  an 
exponential  decay. 

N(t)  =  N0e~,T. 

where  N(t)  is  the  number  of  granule  strings  of  lifetime  r,  and  T 
is  the  mean  granule  lifetime. 

The  granules  found  by  the  center-finding  algorithm  yield  a 
mean  lifetime  of  about  5  minutes,  and  the  granules  found  by 
the  lane-finding  algorithm  have  a  somewhat  shorter  lifetime. 
These  values  are  similar  to  those  of  Kawaguchi  (1980).  who 
followed  granules  from  fragmentation  to  fragmentation. 
Increasing  the  value  of  A r  does  not  increase  the  apparent  life¬ 
time  but  rather  causes  it  to  decrease.  This  is  because  increasing 
the  value  of  A r  causes  the  string  algorithm  to  match  granules 
that  are  close  together  spatially  but  are  otherwise  unrelated. 
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Fio.  1 7  —  Histogram  showing  the  size  distributions  of  granules  found  by  lane  finding  ( solid  lines),  and  by  center  finding  ( dashed  lines) 
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20- — Histograms  of  granule  lifetime*  for  Fourier-filtered  quiet  Sun  data  using  ( top  left )  center  finding  withdr  »  075.  (top  right)  lane  finding  with  dr  =  075. 
i middle  left)  center  finding  with  dr  -  170,  (middle  right)  lane  finding  with  dr  -  170,  (bottom  left I  original  data  using  center  finding  with  dr  »  075.  and  ( bottom  ruthri 
original  data  using  lane  finding  with  dr  -  0"5. 
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Fw.  21.— Mean  size  of  granule*  as  a  function  of  their  measured  lifetime.  The  filled  circles  refer  to  the  mean  of  the  lifetime-average  sue  of  each  granule  during  its 
hfenme.  while  the  filled  squares  are  the  means  of  the  maximum  size  of  the  same  granules. 


For  each  granule  string  we  have  determined  the  lifetime- 
average  and  maximum  size.  Figure  21  shows  a  plot  of  the  mean 
value  of  both  the  lifetime-average  size  (filled  circles)  and 
maximum  size  ( filled  squares)  of  granules  as  a  function  of  their 
lifetime.  It  is  the  result  of  averaging  over  9500  individual 
granule  strings.  The  mean  of  the  lifetime-average  granule  size 
increases  by  nearly  60%  as  the  lifetime  increases  from  30  to  400 
s,  and  then  changes  little  with  lifetime  after  that.  For  the  short¬ 
est  lifetimes  the  means  of  the  lifetime-average  and  maximum 
granule  sizes  are  nearly  the  same.  The  shortest-lived  granules 
have  little  time  to  change  size,  and  may  be  visible  on  as  few  as 
two  images.  After  300  s,  however,  there  is  a  nearly  constant 
difference  of  25%  between  the  mean  value  of  the  lifetime- 
average  size  and  the  mean  value  of  the  maximum  size.  Figure 
22  shows  the  area  of  the  movie  covered  by  features  with  a 
lifetime  f.  In  contrast  to  Figure  18,  it  shows  no  tendency  to 
peak  at  some  characteristic  time. 

We  have  used  the  results  of  our  granule  stride-finding  algo¬ 
rithm  to  study  the  spatial  distribution  of  granules  with  differ¬ 
ent  lifetimes.  Maps  have  been  produced  showing  the  location 
within  our  field  of  view  of  all  granules  with  a  lifetime  within  a 
certain  range.  These  maps  have  been  carefully  compared  with 
the  horizontal  flow  field  by  superposition  with  maps  of  the 
divergence  of  the  flow,  such  as  those  shown  in  Figures  10  and 
11.  In  general  granules  appear  to  be  randomly  distributed  on 
the  surface  of  the  Sun,  independent  of  their  lifetimes.  In  partic¬ 
ular,  the  longest-lived  granules  show  no  preference  for 
occurring  inside  of,  or  on  the  boundaries  of,  mesogranules  and 
supergranutes. 


The  string-finding  algorithm  does  not  yield  granule  lifetimes 
longer  than  those  measured  by  cither  the  original  or  the  sub¬ 
sonic  temporal  ACFs.  Lifetime  results  are  potentially  very  sen¬ 
sitive  to  the  efficiency  at  which  the  string-finding  algorithm 
matches  granules  on  successive  frames.  Given  a  sequence  of 
images  with  10  s  spacing  and  a  granulation  pattern  with  infi¬ 
nite  lifetime,  our  algorithm  would  find  a  lifetime  of  300  s  or  less 
if  its  frame-to-frame  efficiency  were  97%  or  less.  Measuring  the 
efficiency  of  the  algorithm  is  nearly  impossible,  because  of  the 
ambiguity  between  real  misses  and  misses  due  to  granules  just 
being  born,  just  fragmenting,  and  just  dying.  To  get  around  the 
efficiency  problem,  we  have  repeated  the  calculations  allowing 
for  breaks  of  one  and  two  frames  in  the  granule  strings,  so  as  to 
correct  for  noise  in  a  few  frames  that  might  otherwise  artifi¬ 
cially  end  a  granule's  life.  String  mending  did  not  substantially 
modify  the  lifetime  results. 

We  have  also  attempted  to  follow  granules  visually  from 
birth  to  death  on  our  movies.  In  almost  all  cases  examined,  the 
granules  do  begin  and  end  as  the  string-finding  algorithm 
finds.  However,  this  is  not  surprising  after  viewing  the  movies, 
which  show  many  exploding  granules  and  granules  with 
complex  internal  structures.  Below  we  will  discuss  exploding 
granules,  which  may  offer  some  explanation  of  the  apparent 
lack  of  success  of  string  finding,  and  give  some  insight  into  the 
evolutionary  processes  in  the  photospheric  intensity  field. 

Since  our  string-finding  algorithm  can  follow  only  one  of  the 
granule  fragments  when  a  granule  subdivides,  we  have  tried  to 
learn  more  about  granule  lifetimes  by  first  degrading  the 
resolution  to  I'  before  string  finding  in  an  attempt  to  mitigate 
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the  effects  of  fragmentation.  Smearing  was  done  by  appropri¬ 
ate  filtering  in  the  Fourier  domain.  Note  that  the  smeared 
images  are  not  equivalent  to  atmospheric  '‘seeing"  because 
they  are  still  free  of  distortion.  Figure  23  shows  a  histogram  of 
number  of  granules  with  lifetime  t  as  a  function  of  t  for  the 
smeared  images.  Using  the  lane  finder,  the  density  of  granules 
found  is  21  per  100  arcsec1  and  the  mean  lifetime  is  about  10 
minutes.  The  total  number  of  granule  strings  found  is  down  by 
a  factor  of  about  10  when  compared  with  the  original  data. 


Our  string-finding  algorithm  has  been  very  useful  for  mea¬ 
suring  granule  velocities,  because  it  produces  x  and  >- 
coordinates  of  the  center  of  the  granule  throughout  its  lifetime. 
We  have  fitted  a  straight  line  to  each  of  these  sets  of  c  oordi¬ 
nates.  The  slope  of  this  line  is  the  mean  horizontal  velocity  of 
the  granule.  Shown  in  Figure  24  is  the  mean  velocity,  and  the 
standard  deviation  of  the  mean,  as  a  function  of  granule  life¬ 
time.  There  is  an  indication  that  short-lived  granules  move  at  a 
mean  velocity  of  about  1.2  km  s'1.  For  the  longer-lived  gran- 
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tries  the  mean  velocity  is  750  to  850  m  s'*.  This  is  consistent 
with  the  speeds  measured  by  local  correlation  tracking  with 
mask  siam  of  about  ITS,  which  includes  granule  motion  but 
excludes  granule  expansion  and  contraction 

m  EXPLODING  GRANULES 

Granules  that  expand  radially  and  develop  a  dark  central 
region  were  discovered  by  Roach  and  his  coworkers  (Cartier  et 
AL  1908)  and  given  the  name  “exploding  granules."  A  proto* 
typical  example  is  shown  in  Figure  25  (Plate  14).  Usually  the 
expanding  ring  is  not  a  complete  symmetrical  structure  as  seen 
in  the  figure,  but  rather  an  ensemble  of  local  brightenings 
which  are  all  travelling  outward  from  a  center.  They  are  rather 
rare  in  earlier  granulation  data  (MehJtretter  1978;  Namba  ami 
van  Rgsbergen  1977).  but  in  SOUP  data  they  are  common. 
Kitai  and  Kawaguchi  (1979)  have  noted  the  presence  of  a  dark 
central  spot  in  a  large  fraction  of  granules,  but  not  all  of  these 
have  a  redial  expansion  phase.  Referring  to  the  occurrence  of 
exploding  granules,  Beckers  and  Canfield  (1975)  wrote,  "We 
would  not  consider  it  at  aO  unlikely  that  all  granules  would 
show  a  similar  behavior  if  the  quality  of  the  observation  were 
good  enough  "  Recently  Steffen,  Ludwig,  and  Kruss(1988)  cal¬ 
culated  that  the  "exploding" configuration  is  the  steady  state 
form  of  granulation. 

We  have  used  spatial- temporal  maps,  "time-slice"  images, 
to  measure  and  quantify  the  behavior  of  individual  exploding 
granules.  These  images  show  intensity  as  a  function  of  time 
along  a  fine  on  the  Sun.  In  a  time  slice  an  ordinary  granule 
appears  as  a  bright  tube  of  some  length  and  inclination  to  the 


ordinate.  An  exploding  granule  appears  as  a  Y  shape,  and, 
when  the  chosen  line  crosses  the  spatial  center  of  the  exploder, 
the  Y  has  its  maximum  opening  angle.  The  tilt  of  the  arms  of 
the  Y  is  the  radial  expansion  velocity  of  the  exploder.  Measure¬ 
ments  of  the  tilts  show  expansion  velocities  in  the  range  1-2 
km  s*  \  which  is  consistent  with  values  measured  by  Musman 
(1972)  and  Namba  (1986).  Figure  26  (Plate  15)  shows  time 
slices  through  the  center  of  five  exploding  granules. 

We  have  studied  intensively  the  central  32'  x  32'  of  the 
40*  x  40*  quiet  Sun  region,  and  have  found  41  dear  examples 
of  exploders  in  the  1360  s  movie.  The  central  region  was  used 
to  avoid  misidentifying  exploding  granules  occurring  at  the 
edges  of  the  region.  The  average  diameter  of  the  cells  at 
maximum  expansion  is  about  4T2.  Because  it  is  hard  to  recog¬ 
nize  an  event  unless  the  initial  center  or  the  final  ring  can  be 
seen,  a  total  of  about  360  s  at  the  beginning  and  end  of  the 
movie  are  lost  for  detection.  This  suggests  that  at  least  4.0 
exploding  granules  are  born  in  a  10'  x  1(T  area  during  the 
movie,  if  they  were  distributed  uniformly.  This  is  a  lower  limit, 
since,  owing  to  the  relatively  high  density  of  exploders,  some 
events  will  not  be  recognized  because  of  overlap  and  inter¬ 
ference. 

Allen  and  Musman  (1973)  and  Oda  (1984)  have  shown  that 
exploding  granules  are  not  statistically  associated  with  the 
supergranulation  network.  Oda  (1984)  found  some  evidence 
that  growing  granules  may  be  associated  with  roeso- 
granuiation.  Figure  27  (Plate  16)  shows  the  map  of  the  diver¬ 
gence  of  the  horizontal  flow  pattern  overlaid  by  the  location  of 
exploding  granules.  Only  one  exploding  granule  is  well  inside  a 
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region  of  negative  divergence.  The  strong  tendency  of 
exploders  to  occur  in  mcsogranules  at  a  high  rate  is  probably 
the  reason  other  observers  have  noted  that  exploding  granules 
repeat  and  spawn  new  exploders  (Cartier  et  al.  1968;  Kawa¬ 
guchi  1980).  Of  course,  it  also  means  that  all  granules  do  not 
explode. 

Since  regions  of  positive  and  negative  divergence  have 
nearly  equal  areas,  about  eight  exploding  granules  per 
10*  x  10*  area  per  1000  s  are  born  in  cell-center  regions  during 
the  duration  of  the  movie.  The  mean  area  of  maximum  expan¬ 
sion  of  an  exploding  granule  is  about  14  arcsec2,  hence  150 
arcsec2  of  the  average  10*  x  10*  area  are  impacted  by 
exploders  during  the  1360  s  of  the  movie.  This  means  that  in 
every  900  s  or  less  such  regions  are  swept  by  the  expansion 
fronts  of  exploding  granules.  It  is  hard  to  avoid  the  conclusion 
that  virtually  all  the  granules  that  appear  in  the  centers  of 
mesogranulation  cells  either  explode,  are  pushed  aside,  or  are 
terminated  by  an  exploding  granule  in  their  vicinity.  It  also 
means  that  granules  bom  in  the  region  of  negative  divergence 
are  constrained  from  exploding.  In  both  cases,  explosions  are  a 
major  factor  in  determining  the  lifetime  of  granules  in  quiet 
Sun.  These  conclusions  are  verified  by  observations  of  the 
SOUP  and  La  Palma  (Brandt  et  al.  1988)  movies. 

vm.  DISCUSSION 

We  have  attempted  in  this  paper  to  use  computer  algorithms 
and  statistical  techniques  to  identify,  measure,  and  quantify 
properties  of  solar  granulation.  Such  a  plan  is  complicated  by 
the  richness  of  the  phenomena  on  the  solar  surface.  The  inten¬ 
sity  pattern  of  the  solar  photosphere  is  due  to  granulation 
(convective  overshoot  and  turbulent  flows),  oscillations  ( p-  and 
/-modes),  local  internal  gravity  waves,  and  magnetic  fields. 
From  a  single  photograph  it  is  impossible  to  separate  these 
phenomena.  With  the  SOUP  time  sequences  of  images  most  of 
the  effects  of  global  oscillations,  and  some  of  the  effects  of  the 
waves,  can  be  separated  from  the  intensity  pattern.  Simulta¬ 
neous  magnetograms  have  allowed  us  to  isolate  some  of  the 
effects  of  magnetic  fields  as  well.  In  spite  of  the  difficulties  in 
some  areas  the  initial  plan  has  been  successful,  and  we  have 
found  the  following. 

1.  The  temporal  ACF  of  the  photosphere  is  dominated  by 
the  presence  of  solar  oscillations  rather  than  by  the  evolution 
of  granulation.  If  sufficiently  small  areas  of  the  photosphere  are 
studied,  the  ACF  oscillates  with  periods  of  3-6  minutes.  When 
the  solar  oscillations  are  removed  by  Fourier  filtering,  the  tem¬ 
poral  ACF  of  the  photosphere  shows  an  approximately  expo¬ 
nential  decay.  The  decay  lifetime  in  magnetic  field  areas  is 
about  twice  as  long  as  in  quiet  Sun  areas,  which  indicates 
significant  differences  between  the  evolution  of  granulation  in 
quiet  and  magnetic  Sun. 

2.  The  quasi-steady  component  of  the  horizontal  velocities 
of  granules  averages  370  m  s ' 1  (with  4*  mask  size)  in  quiet  Sun. 
The  patterns  in  this  flow  are  on  the  scale  of  meso-  and  super¬ 
granules.  Magnetic  fields  strongly  inhibit  the  average  horizon¬ 
tal  flow  speeds  of  granules,  to  275  m  s' 1  for  weak  fields  and  to 
100  m  s' 1  for  strong  fields.  In  plage  regions  of  nearly  zero 
magnetic  field  the  steady  flows  are  down  by  10%  from  field-free 
quiet  Sun. 

3.  Measurements  of  the  rms  fluctuating  horizontal  velocities 
are  strongly  dependent  on  the  type  of  solar  region  and  the  size 
of  the  aperture  used  for  correlation  tracking  In  quiet  Sun  these 
velocities  go  from  about  0.45  to  1.4  km  s'1  and  in  strong 
magnetic  regions  the  increase  is  from  0.3  to  0.75  km  s' 1  as  the 


measurement  aperture  decreases  from  4’  to  1  ’.  The  velocities 
measured  by  the  smallest  apertures  are  on  the  order  of  veloc¬ 
ities  required  for  models  of  solar  line  widths.  The  rms  velocities 
are  tower  in  the  magnetic  field  areas  by  about  45%. 

4.  The  horizontal  motions  of  granules  alone  are  sufficient  to 
explain  much  of  the  decay  of  the  temporal  ACF.  Thus  tempo¬ 
ral  ACFs  do  not  measure  the  lifetimes  of  granules;  but  rather, 
they  measure  the  evolution  of  the  patterns  in  granulation. 

5.  On  average,  granules  are  brighter  in  magnetic  field  areas 
than  in  quiet  Sun,  and  the  magnitude  of  the  intensity  increase 
is  nearly  independent  of  field  strength  until  pores  form.  Tem¬ 
poral  intensity  fluctuations  are  a  maximum  for  quiet  Sun  and 
decrease  as  magnetic  flux  in  the  local  neighborhood  increases. 
The  fluctuations  decrease  monotonically  with  increasing  field 
strength  until  pores  form.  The  effects  of  magnetic  fields  on 
granulation  extend  well  beyond  (perhaps  2"-4*  or  more)  the 
boundaries  of  the  individual  magnetic  flux  tubes. 

6.  Attempts  to  measure  parameters  such  as  the  mean 
granule  size,  the  granular  density,  and  the  mean  intergranular 
spacing  yield  results  that  are  algorithm-dependent  as  well  as 
resolution-dependent. 

7.  Attempts  to  measure  the  lifetimes  of  all  granules  on  an 
objective  basis  yield  an  average  lifetime  of  about  5  minutes. 
Degrading  the  resolution  of  the  SOUP  images  to  1"  in  an 
attempt  to  reduce  the  effects  of  fragmentation  increases  the 
lifetime  to  about  10  minutes.  The  granules  with  the  shortest 
lifetimes  are  the  smallest  and  have  the  largest  random  horizon¬ 
tal  velocities.  The  longest-lived  granules  appear  to  be  random¬ 
ly  distributed  with  respect  to  the  mcsogranules  and 
supergranules. 

8.  Exploding  granules  tend  to  occur  inside  of  meso- 
granules — in  regions  of  positive  divergence  of  the  horizontal 
flow.  Every  900  s,  on  average,  all  of  the  area  within  these  cells 
has  been  affected  at  least  once  by  the  expansion  fronts  of  an 
exploding  granule.  This  also  means  that  the  evolutionary 
history  of  granules  located  within  the  network  between  meso- 
granules  is  very  different  from  that  of  granules  within  the  cells. 

The  meso  scale  advective  flow  may  be  sufficient  to  account 
for  the  exploders,  or  a  nonuniform  energy  input  associated 
with  the  mcsogranules  may  cause  the  thermal  plumes  in  the 
centers  of  these  cells  to  be  more  vigorous.  Modern  hydrody¬ 
namic  calculations  may  yield  insight  into  the  p.oper  role  of 
mesogranulation. 

In  1962  Leighton  noted  that  there  is  an  asymmetry  between 
the  granules  and  the  intergranular  lanes  visible  on  the  best- 
quality  photographs.  He  then  argued  that  granulation  consists 
of  convective  cells  rather  than  turbulent  flows.  In  1987  Roudier 
and  Muller  determined  that  the  number  of  granules  per  unit 
area  increases  monotonically  with  decreasing  size,  and  they 
argue  that  this  is  evidence  for  a  turbulent  origin  for  smaller 
granules.  Our  observational  results  substantially  agree  with 
both  of  these  authors. 

It  is  easy  to  imagine  that  granulation  results  from  random 
vertical  flows  of  thermal  plumes  penetrating  upward  into  the 
photosphere.  When  the  plumes  reach  the  convectively  stable 
layer,  they  begin  to  expand  horizontally,  become  unstable,  and 
break  up  into  a  turbulent  cascade  with  a  Kolmogorov  spec¬ 
trum.  We  will  resist  forming  conclusions  for  now.  Rather,  we 
will  defer  them  to  further  reports  on  the  detailed  behavior  of 
granules  in  regions  of  convergence  and  divergence  based  on 
ground-based  observations  at  the  new  Swedish  Observatory 
on  La  Palma. 

When  this  work  started,  we  hoped  to  write  the  definitive  life 
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history  of  a  typical  granule.  We  now  know  that  this  is  not 
possible.  The  evolution  of  a  granule  is  strongly  influenced  by 
its  environment,  namely,  by  magnetic  fields  in  its  vicinity  and 
by  its  location  with  respect  to  meso-  and  supergranular  flows. 
Many,  if  not  most,  granules  do  not  evolve  to  completion  but 
rather  are  ended,  probably  randomly,  by  the  radial  expansion 
of  another  granule  in  their  vicinity. 

Small  granules  or  fragments  are  another  problem.  Is  there 
an  objective  way  to  handle  them?  Is  there  a  size  spectrum, 
perhaps  turbulent,  of  small  granules?  Are  the  small  granules 
the  result  of  turbulent  breakup  of  larger  ones?  Clearly,  some 
smaller  fragments  are  generated  during  the  final  stages  of  the 
exploding  granule  process.  Are  other  small  granules  the  signa¬ 
ture  of  a  turbulent  convective  layer?  Are  some  just  small  con¬ 
vective  cells?  Here  again  we  are  tempted  to  make  suggestions, 
but  will  await  the  results  of  new  studies. 

Finally,  we  reemphasize  the  importance  of  the  digital  movie 
techniques  for  studying  and  guiding  data  processing  of  images 
of  granulation  or  any  complex  evolving  system.  Raw  movies 
are  so  chaotic  at  high  resolution  as  to  be  nearly  incomprehen¬ 
sible.  The  three-dimensional  Fourier  filtering  allows  us  to  con¬ 
centrate  on  any  phenomenon  of  interest  which  can  be  isolated 
in  the  (fc,  cu)- plane.  The  all-pervading  nature  of  exploding  gran¬ 
ules  becomes  evident  only  when  the  oscillations  are  removed, 
and  the  differences  between  quiet  and  magnetic  granulation 
also  become  clearer.  The  ability  to  overlay  the  flow  arrows  and 
divergence  maps  on  movie  sequences  was  crucial  to  under¬ 
standing  the  role  of  these  flows.  The  eye-brain  system  is  also 
very  sensitive  to  time  evolution  and  motion,  when  it  is  the 
(datively  narrow  range  the  system  expects,  so  that  it  is  critical 


to  view  movies  at  the  proper  “  brain  speed,"  both  in  forward 
and  reverse  order.  Hopefully,  the  next  few  years  will  see  the 
development  of  image-processing  procedures  that  measure 
more  accurately  what  we  can  see.  Until  that  time,  we  can  only 
comment  that  it  is  indeed  unfortunate  that  journals  cannot 
publish  movies. 
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al  Solar  Observatory  for  providing  magnetograms.  The  SOUP 
instrument  development  and  data  analysis  were  supported  by 
the  National  Aeronautics  and  Space  Administration  under 
contract  NAS8-32805.  The  CCD  camera  used  to  digitize  all 
SOUP  data  was  developed  under  contract  N  ASS-2681 3  for  the 
Coordinated  Instrument  Package  of  the  Orbiting  Solar  Labo¬ 
ratory  (formerly  known  as  the  High  Resolution  Solar 
Observatory).  The  image  processing  systems  for  movie  data 
analysis  using  laser  optical  disks  have  been  supported  by 
Lockheed  Independent  Research  funds. 
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